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We investigate theoretically the nonlinear dynamics induced by an intense pump field in a disor- 
dered planar microcavity. Through a self-consistent theory, we show how the generation of quantum 
optical noise squeezing is affected by the breaking of the in-plane translational invariance and the 
occurrence of photon localization. We find that the generation of single-mode Kerr squeezing for 
the ideal planar case can be prevented by disorder as a result of multimode nonlinear coupling, even 
when the other modes are in the vacuum state. However, the excess noise is a non-monotonous 
function of the disorder amplitude. In the strong localization limit, we show that the system be- 
comes protected with respect to this fundamental coupling mechanism and that the ideal quadrature 
squeezing generation can be obtained. 

PACS numbers: 42.50.Dv, 71.55.Jv, 42.65.Hw, 42.50.Pq 



Recently, the intrinsic properties of optically active dis- 
ordered media have been attracting a considerable inter- 
est. In particular, an impressive deal of research has been 
devoted to the study of the interplay between lasing and 
photon localization in random media (THE!, described in 
the framework of semiclassical Maxwell- B loch equations. 
Concerning the impact of disorder on purely quantum op- 
tical properties, research appears to be in its infancy. Re- 
cent works have investigated how an input squeezed light 
propagate through a linear random medium [§[ or how a 
slight nonorthogonality of cavity eigenmodes affects the 
squeezing Q . The generation of quantum optical squeez- 
ing for quantum applications has been stimulating both 
in atomic and condensed matter physics @,@|. In the case 
of an active medium with a third-order optical nonlinear- 
ity, it is known that the coherent resonant pumping of a 
photonic mode can produce a Kerr quantum squeezing 
of the optical field [Io|, [HI • To understand mechanisms 
preventing the ideal squeezing generation is clearly a fun- 
damental issue. To the best of our knowledge, the pump- 
induced generation of quantum squeezing in a nonlinear 
disordered medium has not been explored. 

Here, we wish to explore the role of photon localization 
on the quantum squeezing generation of a 2D nonlinear 
photonic system (like a planar microcavity) in presence 
of a disordered potential. We have consistently deter- 
mined the nonlinear quantum dynamics of exact photonic 
modes in presence of disorder. Our results shows that 
the multimode nonlinear coupling is responsible for an 
excess noise, which can prevent the quantum squeezing 
generation. Interestingly we have found that the excess 
noise depends non-monotonously on the disorder poten- 
tial amplitude AV. The maximum of noise is obtained 
for AV comparable to the homogeneous linewidth of the 
mode. However, increasing further the disorder, we find 
that the ideal quantum squeezing generation can be even- 
tually recovered, as the pumped photon mode becomes 



"protected" from the multimode coupling in the strong 
localization limit. 

Let us consider a model quantum Hamiltonian for a 2D 
photonic system (like a planar microcavity) subject to a 
third-order nonlinearity and to an in-plane (disordered) 
potential for the the photonic motion : 
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where uok Z) k = ^ yfef + A? is the frequency of the pho- 
tonic mode (fe is the in-plane wavevector, k z is the quan- 
tized vertical wavevector, while e is the dielectric con- 
stant). Vk-k' is the Fourier transform of the in-plane 
(disordered) spatial potential affecting the in-plane mo- 
tion of photons, while g is the nonlinear coupling coeffi- 
cient for the third-order optical nonlinearity. Moreover, 
dk is the annihilation operator of a cavity photon with 
in-plane wavevector k. 

If we diagonalize exactly the Hamiltonian part includ- 
ing the free cavity photon motion and the disordered po- 
tential, we can rewrite the Hamiltonian as: 
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where uji and di are the eigenfrequency and annihilation 
boson operator of the z-th photon eigenmode in pres- 
ence of the disorder potential k,l G {0,1,2,3,...}). 
The quantity gijki is the nonlinear coupling coefficient 
depending on the eigenmode normalized wavefunctions 
i/ji(r). Namely, it reads: 
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A continuous wave coherent optical pump applied to 
the system is able to drive coherently the intracavity pho- 
tonic modes and produce a net mean- field (hi). In the 
Heisenberg picture, the quantum statistical properties 
of the excited light depend on the fluctuation operator 
8ai = ai — (ai) . By linearizing with respect to the fluctu- 
ations and by considering the coupling to the extracavity 
field and to non-radiative reservoirs, we get: 
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and then Eq. (j3J) is reduced to 
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The dynamics of the intracavity fluctuation operator de- 
pends on the mode eigenfrequency U0i and on the ra- 
diative loss rate r* ad and the non-radiative one r^ R . 
The terms proportional to gij^i describe the nonlinear 
renormalizations and the multimode coupling. The op- 
erator Fi is the radiative Langevin force for the i-th 
mode and SFi = Fi — (Fi), where (Fi) is the excita- 
tion amplitude induced by the coherent optical pump- 
ing, where (SFi(t)SFj (£')) = SijS(t - t')r™ d . The non- 
radiative Langevin force is such that (F^ R (t)F^ R ^ (t f )) = 
SijS(t - t , )r i NR (Nl h + 1), where N* h is the thermal pop- 
ulation of i-th mode. 

The quadrature amplitude of the i-th intracavity mode 
is defined as Xfy(t) = ai(t)e [uj ^- w + a\(t)e-' l ^ t+ie . The 
quantum noise is then quantified by the correlation func- 
tionCft(r) = (X™(t)X%(0)) - (X$(t))(X%(0)) = 
(5Xl n e (r)5Xl n e (0)). The intracavity quadrature noise is 
given by the equal-time correlation function CYq{t = 
0). The extracavity quadrature is defined as X^(t) = 
& ouT^ p t-ie + & t ? UT t(j) e -i"pt+tf where af VT (t) is 
the output operator with the same transverse profile of 
the i-th intracavity mode (see Supplementary Materi- 
als). The extracavity quantum noise is quantified by 
the Fourier transform of correlation function Sffffa) = 
/~ dr e iwT (r )Xff (0) ) - (Xff(r))(Xff(0))}. For 

a normal vacuum or for a coherent state, C-^fY = 0) = 
Sffffa = 0) = 1 (standard quantum limit). The in- 
tracavity and extracavity squeezing are obtained when 
C i n e( r = 0) < 1 and S°ft(oj = 0) < 1, respectively. 

In the case of a perfectly planar cavity (no disorder), 
the eigenfunctions are planewaves ^j(r) = e lkj V /\/~S and 
the nonlinear coefficients are giju = gfiki+kj,k k +ki • If 
we consider the optically pumping of the cavity by a res- 
onant laser with in-plane wavevector q = 0, only the 
ground intracavity mode (feo = 0) is coherently excited, 
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Therefore, since Eq. (j5]) has no interaction with the fluc- 
tuations of the other modes, the dynamics reduce to 
that of the well-known single-mode Kerr model. In the 
case of negligible non-radiative damping and no thermal 
population (A^ th = 0), the minimum quantum noise is 
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obtained for pumping intensities close to the bistability 
turning point [12l4l6|. It is known that such extracavity 
perfect squeezing is due to an interference effect between 
the incident pump field and the reemitted light. 

In presence of disorder, the situation is qualitatively 
different. We have solved Eq. (j4]) by considering a ran- 
dom potential with amplitude Ay and correlation length 
A c . We have calculated the squeezing of the ground intra- 
cavity mode (denoted by i = 0) assuming that only the 
ground mode is coherently excited ((Fi) = Si^FQ Ump (t)). 
This means that only the ground mode has a coher- 
ent amplitude (clq) and that the expectation values in 
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Eq. dH) are g 

and (ctidj) = (^0^,0 (&o) 2 • Following Eq. (j4|, the ground 
mode (i = 0) is coupled to the higher modes via the terms 

2#oooj(^o^o)^j and #ojOo(^o^o)<^J- Such intermode cou- 
pling terms produce an excess noise for the ground pho- 
tonic mode. 

Representative results are depicted in Fig. [TJ We have 
considered parameters for a GaAs microcavity (details in 
the caption). For the case N\ h = (no thermal popula- 
tion) , the only relevant damping parameter for the intra- 
cavity dynamics is the total loss rate r[ ad + r^ R = r. 
Of course, for the extracavity noise statistics, the ratio 
rQ R /r™ d is crucial. In the calculation of Fig. [TJ we 
have considered rQ ad ^> r^ R (ideal condition). In the 
case A^ h > 0, the non-radiative loss rate plays an impor- 
tant role. For the modes i ^ 0, we have considered for 
simplicity r^ R = r. We have considered a finite-size 2D 
system (100 /imx 100 jim) with periodic boundary condi- 
tions and the correlation length of disorder is A c = 5 //m. 
The pump frequency is uj p = uJi=o + F. In Figs. H{a) 
and[]Jb), we present results for the minimum intracavity 
quadrature noise C]= and extracavity one S°=q, respec- 
tively, as a function of the disorder potential amplitude 
AV, normalized to the homogeneous broadening r. The 
solid line refers to the case where N^ = 0, i.e., the non- 
pumped modes are in the vacuum state with no excess 
thermal noise. For negligible disorder (AV/T <C 1), the 
minimum intracavity noise takes the value C£L = 0.5, 
and the extracavity noise is S?=q = 0. For AV/T = 0.1 
(point (i)) the intracavity and extracavity noises are only 
slightly increased to 0.55 and 0.13, respectively. Increas- 
ing AV/r, the minimum noises significantly increase and 
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FIG. 1: (a) Minimum value C\= of the quadrature quan- 
tum optical noise of the ground intracavity mode versus the 
disorder potential amplitude AV (normalized to the mode 
homogeneous broadening r). The dotted line (Cl= = 1) 
is the standard quantum limit. The solid line is calculated 
under the condition that only the ground mode is driven co- 
herently, while all the other modes are in the vacuum state. 
The dashed line is obtained by considering a thermal popu- 
lation iV^o = 1, while only the ground mode is coherently 
pumped, (b) Minimum value >S°=o of the quadrature noise 
of the extracavity field which has the same transverse pro- 
file as the ground intracavity mode, (c) Surface plots of the 
disordered potential (bottom) and of the ground mode pho- 
tonic wavefunctions (top) for the different amplitudes (i), (ii), 
(iii) of the disorder potential, also indicated in (a) and (b). 
Parameters: 100/mi x 100/mi system with periodic bound- 
ary conditions; disorder potential with correlation length 
A c = 5 jLim; huk z ,o =_ 1.5 eV; e = 13; g = 0.01 meV/xm 2 ; 

= 0; 



nr 



0.1 meV; hr?$ 



+ r. Averaging over 100 disorder configurations. 



reach a maximum values of C£L = 0.77 and S£=q = 0.75 
for AV/T = 2 (point (ii)). Increasing AV/T further, the 
noises decrease and recover the values Cf^ Q = 0.5 and 



s° 



(for example, AV/T = 100, point (iii)). In 



Fig. H{c), we show a surface representation of the pho- 
tonic ground mode squared wavefunction |^o(r)| 2 (top) 
and of the disordered potential (bottom), showing the 
increasing degree of localization of the photonic ground 
wavefunction with increasing value of A V. Note that the 
dashed line in Figs. QJa) andQJb) instead is for the case 
with N^ = 1, meaning that the non-pumped modes 
have excess noise due to a thermal population. In such 




Disorder amplitude AV / r 

FIG. 2: (a) Localization length Lo of the ground mode nor- 
malized to the system size L p i ane versus the normalized dis- 
order amplitude AV/T. (b) Solid line: effective number 
iVeff = 7rrVKDOs(cJo) — 1 of modes weighted by the nonlin- 
earity to the ground mode. Dashed line: actual number 
^Vact = 7ri^pDOs(<^o) — 1 of modes, by taking into account 
the normal density of states. 



a case, there is a similar behavior to the case N^ Q = 
(solid line), but, importantly, the minimum noise can be 
larger than 1 (larger than the standard quantum limit). 

In order to understand the origin of this interesting 
behavior of the nonlinear quantum squeezing generation 
in presence of disorder, we have studied in detail the 
localization properties, the density of states of the dis- 
ordered eigenstates and the nonlinear coupling between 
them. We have found that the generation of quantum 
squeezing can be forbidden if (a) there is a significant 
number of modes close in energy to the ground mode 
{uJi>\ — uoq ~ r) and (b) they are significantly cou- 
pled to the ground mode (goooi has to be comparable 
to goooo)- In order to quantify the localization, we have 
considered the localization length Lo = y (\r — (r}\ 2 ) of 
the photonic ground mode, where the average (• • •) is 
taken on the ground mode state. In Fig. Efa), we plot 
Lo divided by L p \ ane (the size of a planewave mode) as 
a function of disorder amplitude AV. For small values 
of AV, Lo/Lpiane — I an d, since the translational in- 
variance is broken only very weakly, the nonlinear cou- 
pling coefficients are such that goooi/goooo — ^,o- How- 
ever, increasing A V, the wavefunction becomes more and 
more localized and a transition from Lo/L p \ ane = 1 to 
Lo/Lpiane ^ I occurs. In the strong localization limit, 
we have found that an energy gap between the ground 
photonic mode and the higher modes is opened (see the 
sketches of energy structures in Fig. Efa)). This en- 
ergy gap effectively decouples the ground mode (i = 0) 
from the other modes (i > 1) and protects the quan- 
tum squeezing generation from the multimode coupling. 



Hence, the ideal values for single- mode Kerr squeezing 
are recovered both in the case of negligible and very 
strong localization, albeit for very different reasons. In 
the case of no disorder, the selection rule due to the 
translational invariance is responsible for the mode de- 
coupling; in the strong localization limit, the decoupling 
is due to a frequency gap 0Ji>i — ujo ^> T. In between, 
the multimode coupling produces an excess noise, which 
is maximum for a critical value of AV/T (around 2 for 
the considered type of disorder). 

In order to get further insight into the nonlinear mul- 
timode coupling, we have considered the density of pho- 
tonic states Pdos(^) = J2i where £i(u) is a nor- 
malized Lorentzian function (J duo Ci(uo) = 1) with cen- 
ter uji and width iy ad + r^ R . Having the density of 
states Pdos(<^), it is possible to compute the number 
N ac t = 7r J TpDOs(^o) — 1 of modes whose eigenfrequencies 
differ from ujo less than r. Such a number 7V ac t is plotted 
in Fig.[2fb) (dashed line). It is apparent that 7V act is not 
responsible alone for the non-monotonous dependence of 
the quantum noise versus disorder amplitude in Fig.QJa). 
Indeed, as apparent from Eq. (|4]), the Kerr quantum 
squeezing is generated by the term #0000(^0^0)^0? while 
the excess noise due to the coupling to the i-th mode 
(i 7^ 0) comes from the term goioo(a>oa>o)5al • Hence, the 
ratio between the amplitudes goooi and #0000 plays an 
important role. It is therefore instructive to study the 
behavior of the quantity 
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^dos(^) is a modified density of states where the con- 
tribution of each state is weighted by its nonlinear cou- 
pling to the ground mode. Such simple weighted den- 
sity of states indeed qualitatively reproduces the non- 
monotonous dependence of the excess noise as a func- 
tion of the disorder amplitude in Fig. []Ja). This is 
confirmed by looking at the effective number of modes 
obtained with the weighted density of states, namely 
Af e ff = tt^Wdos^o) — 1? which is plotted by a solid line 
in Fig.E{b). 

From the above discussion, it is now clear that, in pla- 
nar cavities with disorder amplitude AV comparable to 
the mode homogeneous broadening r, the generation of 
quantum squeezing can be severely limited by nonlinear 
coupling between disordered states. This critical condi- 
tion A V c± r is actually achieved in semiconductor GaAs 
microcavities, which have been experimentally investi- 
gated in recent works HI, [13, [Hj. However, our results 
also show that the system where photonic energy levels 
are separated much more than the homogeneous broad- 
ening r the ideal single-mode Kerr squeezing can be re- 
covered. Therefore, in this context, it appears natural 
to consider artificially localized photonic systems such as 
micropillars [l8| to achieve a protected quantum squeez- 
ing generation. In GaAs-based microcavities, discrete 



photon levels are obtained already in a micropillar with 
a 4 /im diameter and a level separation larger than the 
broadening can be achieved [l9j. Of course, the complete 
suppression of extracavity noise can be detected only by 
using a local oscillator [lOj with the same spatial shape 
as the ground photonic wavefunction. 

In conclusion, we have presented a theory for the role of 
photon localization on the nonlinear generation of quan- 
tum squeezing in a disordered planar microcavity. We 
have shown that the ideal single-mode quantum noise 
reduction is limited by nonlinear multimode coupling be- 
tween disordered eigenmodes, even when the modes are 
in the vacuum state. Protection from multimode cou- 
pling is obtained in the strong localization limit [20] . Our 
study shows the interesting link between nonlinear quan- 
tum noise properties and photon localization and may 
stimulate a wide range of experimental and theoretical 
investigations in a variety of systems. For example, the 
link between Anderson localization and quantum noise 
in 3D photonic systems Q, or the quantum noise statis- 
tics of superfluid light in a nonlinear disordered medium 
22 . [23| are intriguing problems to explore in the future. 

We wish to thank A. Amo, A. Bramati, I. Carusotto, C. 
Fabre, E. Giacobino, J. Ph. Karr for discussions and/or 
for a critical reading of the manuscript. 



SUPPLEMENTARY MATERIAL 
1. Calculation of in-plane eigenmodes 

As disordered potential Vh for the in-plane motion of 
photons, we have considered a Gaussian function 
with amplitude AV and correlation length A c . Namely, 
we have taken Vh = AVe~ Xc k / 4 + l6>k , where the phase 
#k is given by a random function satisfying Oh = —0-h to 
make V(r) = ^2h Vk£ lk r real. The 2D system with size 
L x x L y is divided into a N x x N y grid (N x = N y = 128 in 
our numerical calculations) and the in-plane wavevector 
k = (k Xl k y ) in the first Brillouin zone is used: k XyV = 
2irn/L x , y and n = -N x , y /2, -N x , y /2 + 1, . . . , N x , y /2- 1. 
The eigenfrequency uoi and function ^i( r ) m the 2D 
system is calculated by solving a eigenvalue problem 
[huj kz (-\V) + V(r)]^i(r) = hjJitj)i{r), where uj kz (-\V) = 



V k z 2 — V 2 . We have used periodic boundary condi- 
tions for the photonic wavefunction. In the numerical 
diagonalization, the matrix is relatively sparse. 



2. Treatment of extracavity fields 

For the treatment of extracavity fields in Eq. (4) on the 
main paper, we have followed a generalized input-output 
approach for the present system. The coupling between 
the intracavity and the extracavity fields is described by 
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the following Hamiltonian: 



ifi 



[ K i,q<*i<Xq(u) ^ q OL q {uj)ai 



(7) 

a q (uj) is the annihilation operator of the extracavity field 
with frequency uj and in-plane wavevector q and the cou- 
pling coefficient is represented by the Fourier trans- 
form ijji^ of eigenfunction ipi(r) as = V r Tad ?pik. 
Here, it is convenient to consider a set of extracavity 
fields with the same transverse profile as the intracavity 
eigenmodes, namely: 



In terms of this basis, Eq. © is rewritten as 
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Therefore, the extracavity mode with the same trans- 
verse spatial profile as the i-th intracavity mode is not 
affected by the other orthogonal intracavity modes, and 
then Eq. (4) on the main paper and the input-output 
relation 



;OUT 



(10) 



are derived by the standard approach |8|, |9j. a| N ^ OUT (co>) 
is the so-called input /output operator (operator &i(u), to) 
where to — >• — oo/oo) is the initial/final time of the dy- 
namics, and the Langevin force operator in Eq. (4) on 
the main paper is written as 



Fi(t) = (Fi(t)) + SF(t) 



— f 



(11) 

In practice, the dynamics of extracavity field &i(cj) can 
be measured by using a local oscillator spot with the ap- 
propriate shape, if the wavefunction is mostly con- 



centrated within the light cone. Actually, in our calcu- 
lation, there is negligible amplitude of the fundamental 
mode ipo,q outside the light cone. 



3. Details about correlation functions 

The quadrature amplitude of the z-th intracavity mode 
is defined as X^ e (t) = ai(t)e [uj ^- [G + a\(t)e-' luJ * t+ie . The 
quantum noise is then quantified by the correlation func- 
tion C$(r) = (X™(t)X§(0)} - (X$(t))(X§(0)) = 
(6X$(t)5X§(0)), where 5X§ = Xfy - (Xfy). The 
intracavity quadrature noise is given by the equal-time 
correlation function C^(r = 0). The extracavity 
quadrature is defined as X?f(t) = a? VT (t^^ + 
a? UT ^(£)e _lu;p * +10 . Its quantum noise is quantified by 
the Fourier transform of correlation function S^ l (u}) = 
/~dr e^[(Xff(r)Xff(0)) - (X°f(r))iX°£(0))]. As 
already well explained in the literature 
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the ex- 
tracavity field quantum noise properties are affected by 
the interference between the incident (input) field and 
the reemitted light by the cavity. If the quadrature noise 
of the fundamental intracavity mode is Q= = 0-5? it 
can be shown that a perfect squeezing (S£=q = 0) can be 
obtained outside in the case <C rQ ad . 



4. Linearization Bogoliubov method 

Concerning the stability of the mean-field solutions 
(cbi), we note that the brackets in Eq. (4) on the main 
paper are calculated as (a\aj) = ^,o^j,o|(^o)| 2 + ^i,j^i h 
and {(liCLj) = #i,o<^,o(^o) 2 - Moreover, we treat |(ao)| 2 as a 
variable proportional to the pump intensity and check the 
system stability by a linearization Bogoliubov method. 
When we consider TV eigenmodes of in-plane photonic 
motion, the size of the Bogoliubov matrix L is 2N x 27V. 
When L is divided into N x N blocks with each size 2x2, 
the block is written as 



/ 7-rrad , j-iNR 

({Soj - i £i + r > )5 ij + 2g l00j \(ao}\ 2 
V -9t jo0 {ao)* 2 



j-i r a d _i_ j-i N FL 

(Su - i « 2 * - 2 5 * 00i |(a )| 2 / 



(12) 



where Sou = uji + 2giooiN^ h — uj p . The system is sta- ing the pump intensity in the stability region. For the 

ble when the imaginary part of its eigenvalues are all system with only one photon mode, the minimum val- 

negative. In order to find the minimum of intracavity ues C*= = 0.5 and = are obtained at two pump 

and extracavity quantum noise, we calculated the val- intensities at an instability boundary. However, for the 

ues C-" = min[C^ e (0)] and S£=q = min[5°^Q e (0)] system with many modes and no translational invariance, 

obtained by tuning the quadrature phase 6 and by tun- in general, the system has many instability regions and 
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the minimum is obtained not just at the boundaries. We 
considered the minimum values below the first instability. 
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